The height evolution of the 'true' CME mass derived from 
STEREO COR1 and COR2 observations 

B.M. Bein 1 , M. Temmer 1 
1 Kanzelh6he Observatory-IGAM, Institute of Physics, University of Graz, Universitatsplatz 

5, 8010 Graz, Austria 

A. Vourlidas 2 

2 Space Science Division, Naval Research Laboratory, Washington, DC, USA 

A.M. Veronig 1 

1 Kanzelh6he Observatory-IGAM, Institute of Physics, University of Graz, Universitatsplatz 

5, 8010 Graz, Austria 

and 

D. Utz 3 - 1 

1 Kanzelh6he Observatory-IGAM, Institute of Physics, University of Graz, Universitatsplatz 

5, 8010 Graz, Austria 

3 IAA - Instituto de Astrofisica de Andalucia, CSIC, Glorieta de la Astronomfa, s/n, 18080 

Granada, Spain 

Received / Accepted 



Received 



accepted 



-2 - 
ABSTRACT 

Using combined STEREO-A and STEREO-B EUVI, COR1 and COR2 data, 
we derive deprojected CME kinematics and CME 'true' mass evolutions for a 
sample of 25 events that occurred during December 2007 to April 2011. We 
develop a fitting function to describe the CME mass evolution with height. The 
function considers both the effect of the coronagraph occulter, at the beginning of 
the CME evolution, and an actual mass increase. The latter becomes important 
at about 10 to 15 R and is assumed to mostly contribute up to 20 R . The mass 
increase ranges from 2 to 6% per R and, is positively correlated to the total 
CME mass. Due to the combination of COR1 and COR2 mass measurements, we 
are able to estimate the 'true' mass value for very low coronal heights (< 3 R ). 
Based on the deprojected CME kinematics and initial ejected masses, we derive 
the kinetic energies and propelling forces acting on the CME in the low corona 
(< 3 R ). The derived CME kinetic energies range between 1.0 — 66 • 10 23 J, and 
the forces range between 2.2 — 510 • 10 14 N. 

Subject headings: Sun: coronal mass ejections, Sun: activity 
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1. Introduction 

The Sun's atmosphere is losing mass continuously via the solar wind outflow and 
sporadically via coronal mass ejections (CMEs). In this paper we focus on the mass ejected 
via CMEs. CME velocities derived from coronagraphic observations range from ~100 
to 3000 km s- 1 



le.g. 



Yashiro et al. 2004; Gopalswamy et al. 2009). It was found that 



the CME peak acceleration may occur very low in the corona (< 1.5 R in heliocentric 
distances; e.g. Maricic et al.||2007 ; Bein et aJ7||2011 ) . It is therefore necessary to consider also 
non-coronagraphic observations for the study of the CME impulsive phase. We combine 
EUV images with coronagraphic observations to obtain detailed CME kinematics from its 
initiation close to the solar surface up to about 15 R , including all three phases of CME 



evolution (initiation, impulsive acceleration, propagation phase) described in Zhang et al. 



(2001, 2004). 



Vourlidas et aT| ( pOlO) analysed a sample of 7668 CMEs observed by the LASCO/SOHO 
coronagraphs during solar cycle 23 and found that the distribution of CME masses ranges 
from 10 13 to 10 16 g, with a few outliers up to 10 17 g and as small as 10 11 g. When observing 
CMEs from one spacecraft, i.e. from only one vantage point, we observe the projection of 
a three dimensional structure on the plane-of-sky (POS) and the derived CME quantities, 
like the CME velocity, acceleration or mass are underestimated. The larger the deviation of 
the propagation direction of the CME from the POS, the larger the projection effects. With 
the advent of the STEREO mission it is now possible to derive propagation directions and 



deprojected CME quantities using various methods (e.g. Mierla et al. 2010). Triangulation 



methods estimate the direction of CMEs by combining height-time curves derived from two 



vantage points (e.g. de Koning et al. 2009 Liewer et al. 2009 Temmer et al. 2009). In the 



forward modeling method applied to stereoscopic observations, a flux rope is fitted to the 
white light observations, which enables us to estimate the CME propagation direction as 



well as its 3D geometry (Thernisien et al. 2009). Colaninno & Vourlidas (2009) presented a 
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technique to estimate the 3D, i.e. deprojected, CME mass, by assuming that the difference 
in the observed mass from two vantage points is solely due to projection effects. We would 
like to note that although stereoscopic observations are used, we still do not know the real 
3D structure of the CME body, and therefore, the 'true' mass remains uncertain. Similar 
to the triangulation method, the masses observed from both vantage points are corrected 
for different propagation directions until both mass estimates yield the same result. Thus 
also the CME propagation direction can be derived with this method. 
Based on these techniques we present a study of deprojected CME quantities for a set of 
25 events observed by both STEREO spacecraft. The set of CMEs was chosen according 
to observational criteria: CMEs must be well observed in both STEREO spacecraft with 
the instruments EUVI, COR1, and COR2. With this we are able to study the CMEs 
starting from the solar surface up to 15 R Q and to derive detailed deprojected height-time, 
velocity-time and acceleration-time curves. We show how COR1 and COR2 measurements 
can be combined to track the CME mass evolution over a large height range. Due to the 
occulter disk of coronagraphs, measurements of the early evolution of the CME mass do not 
reflect the real CME mass accurately. Therefore we develop a fit function which considers 
both geometric and physical changes in the CME mass evolution. From this fit function 
we derive the mass at the earliest observational evolution of the CME as well as define a 
final CME mass, i.e. the CME reaches a constant mass above some height. Furthermore we 
present deprojected CME kinetic energies and forces in the low corona. 



2. Data 



The Sun-Earth-Connection Coronal and Heliospheric Investigation (SECCHI) package 
(Howard et al. 2008) onboard the Solar Terrestrial Relations Observatory (STEREO) 
mission (Kaiser et al. 2008) provides stereoscopic CME observations from the Sun up to 
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1 AU, by the combination of five different instruments: one Extreme Ultraviolet (EUV) 



Imager (EUVI; Wuelser et al. 2004), two white light coronagraphs (COR1 and COR2) and 
two white light heliospheric imagers (HI1 and HI2; |Eyles et al.|2009j ). For the present study, 
CORf and COR2 coronagraphic observations in combination with EUVI data are used to 
derive 3D CME kinematics and mass evolutions. 

With a field-of-view (FoV) of 1.7 R & , EUVI observes the upper chromosphere and lower 
corona in four different wavelength bands. For our study we use 195 A images with a time 
cadence of 2.5 to 10 minutes. COR1 observes the inner corona between 1.4 to 4 R Q with 
a time cadence of 5 to 10 minutes. The FoV of COR2 observations ranges from 2.5 to 15 
Rq. We use total brightness images either obtained directly or derived from polarization 
sequences to achieve the maximum COR2 cadence of 15 mins. For 3 events in our sample 
we have lower time cadences of 30 to 60 minutes. The analysis comprises a set of 25 CME 
events which occurred between December 2007 and April 2011. 



3. Analysis 

For the sake of brevity, we demonstrate our data analysis procedures using a single 
example, the 2010 April 3 CME. 



3.1. Kinematics 



We measure the height-time profile of the CME leading edge in EUVI, COR1 and 
COR2 images to derive detailed CME kinematics starting from the solar surface up to 
about 15 Rq using the methods described in Bein et al. (2011). Figure [I] shows an image 
sequence for the 2010 April 3 event, observed by EUVI, COR1 and COR2 from STEREO- A 
and STEREO-B. In each image we determine the leading edge (black solid lines in Fig. 
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[T| and the plane-of-sky (POS) propagation direction (black radial lines in Fig. [TJ, which 
cross at the center of the CME front. The projected CME height is measured from 
STEREO-A as well as from STEREO-B observations. As reference point we identify from 
flaring signatures in EUVI data the location of the source region of the CME at the solar 
surface. From this reference point we calculate the average height for all measurement 
points on the CME leading edge, that are within ±5° of the POS propagation direction. 
By applying the forward modeling of a flux rope structure to STEREO-A and STEREO-B 
coronagraphic observations, we obtain the propagation direction of the CME (Thernisien et 
al., 2009). This direction is used to correct the height-time curve in order to get deprojected 
kinematics. 

Figure [2] shows the resulting deprojected kinematics for the 2010 April 3 event. The top 
panel shows the height-time curve, the middle panel the CME velocity and the bottom panel 
the acceleration profile. Our typical measurement error is 0.03 R & for EUVI, 0.125 R & for 
COR1 and 0.3 R Q for COR2. Our best estimate of the uncertainties due to the measurement 
in different wavelengths (EUV, white light) across the instruments are included in the 
measurement errors. From a spline fit on the derived height-time curve, we obtain smooth 
velocity and acceleration profiles from the first and second derivative, respectively. The 
maximum value in the CME acceleration profile is used as peak acceleration a max . The 
peak velocity v max is defined as that value in the velocity profile, where the acceleration 
has decreased to 10% of its maximum value. The grey shaded area in the velocity-time and 



acceleration-time plot represents the error range on the spline fit. For more details see Bein 
elldlpmj). 
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3.2. Mass measurements 



Applying the technique described in Vourlidas et al. (2010) on COR1 and COR2 



observations, we derive the CME mass evolutions from both STEREO vantage points. Base 
difference images are constructed by subtracting a pre-event image, containing no CME 
signatures or other disturbances. Both the pre-event image and the CME image itself are 
corrected for instrumental effects and calibrated in units of mean solar brightness (e.g. 



Poland et al. 1981). In the ideal case this procedure results in an image, which shows a 
brightness excess caused by the CME only. The brightness excess is converted into electron 
excess with the Thomson scattering formulation assuming that all electrons are located 
on the spacecraft POS and a composition of 90% H and 10% He. Because the Thomson 
scattering along the line-of-sight (LOS) is the strongest for electrons located on the POS, 



this method leads to an underestimation of the mass of the order of a factor of 2 (Vourlidas 



et al. 2010). Depending on the width and the propagation direction of the CME, this can 
affect the entire CME or parts of it. 

The second to fifth row in Fig. [T] show examples of mass images, in which the CME regions 
are outlined by using the 'sector' method. We define in each image a sector, which contains 
the CME structure. For all observations (COR1 and COR2) from the same STEREO 
vantage point and relating to the same event, the width of the CME sector is kept the 
same. For COR1 the lower boundary of the sector is set at 0.1 R above the occulter, for 
COR2 at 0.2 R above the occulter to avoid the increased noise around the occulter disk. 
The upper boundary of the sector is dependent on the height of the CME leading edge. 
Because each pixel value of the mass image indicates the mass along the LOS at that point, 
we obtain the CME mass by summing up all the pixel values within the defined sector. 
Assuming that the same mass is observed from both spacecraft, mass measurements from 
STEREO-A and STEREO-B are corrected for different CME propagation directions until 



a best match is found (Colaninno & Vourlidas 2009). This method is applied on the entire 
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time series of images to derive a stable propagation direction. A fixed value (mean value of 
the propagation direction derived from the last three C0R2 observations) is used to derive 
the 3D mass evolution for COR1 and COR2 observations. 



3.3. Error Analysis 

Several factors may affect our mass measurements. To get a lower and upper limit of 
the CME mass measurements, we compare for several sample events the results derived 
from the region of interest (ROI) method with the results derived from the 'sector' method. 
The ROI is manually selected enclosing the boundary of the CME structure. In cases, 
where the boundary of the CME is blurred, we defined the ROI very tight, to get a lower 
limit of the CME mass measurements. The 'sector' method includes all parts of the CME 
but sometimes also coronal disturbances (e.g. streamer deflections), so it serves as an upper 
limit. The differences in the absolute mass estimates derived from the two methods are 
on average 30%. The error due to subtracting different pre-event images, which differ 
in time up to 13 hours, is found to be of the order of 10-20%. Generally we select our 
pre-event images very carefully, thus the actual uncertainties are much lower. To consider 
the maximal possible error we apply for all events error bars of ±15% of the current mass 
value on our measurements. 



Vourlidas ct al. (2010) discussed several other effects influencing mass measurements in 
LASCO images (instrumental effects, coronal background, composition of the coronal 
material) and showed that they are in the order of some percent, i.e. within our error 
estimate. In contrast, projection effects can lead to an underestimation of the mass of a 
factor of two. 

We minimize these projection effects with the calculation of the 3D mass evolution. 
However, there are still some uncertainties in the determined propagation direction. The 
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propagation direction, which is used for calculating the 3D mass, was taken from the mass 
calculation method. To estimate an appropriate error, we also use propagation directions 
derived from other methods (triangulation, forward modeling) and compare the masses 
derived with the different directions. Only for five events out of our study these differences 
are larger than 15% (~20% for 2008 April 26, 2010 August 2; ~30% for 2010 September 2, 
2011 January 31, 2011 March 21). Compared to the rest of the sample under study, those 
events were found to be either very wide and faint (2010 September 2, 2011 January 31) 
or that their propagation directions deviate more than 50° from the spacecraft POS (2008 
April 26, 2010 August 27, 2011 March 21). 



4. Results 

Figure [3] shows COR1 and COR2 mass measurements derived from STEREO-A (top 
panel) and STEREO-B observations (middle panel) and the 3D CME mass evolution 
(bottom) against time. COR1 and COR2 observations at the same time do not match 
because they do not include the same CME volume due to the different occulter sizes of 
the instruments. For example, a CME observed in COR2 below ~5 R shows less mass 
than observed in COR1 for the same height range. This occulter effect has been noted since 



the very first LASCO mass measurements (Vourlidas et al. 2000). To test and properly 



quantify the occulter effect, we calculate the COR1 mass applying the COR2 occulter size 
and position (cf. Fig. [4]). These measurements are plotted as green triangles in Fig. [3] and 
fit well with the actual COR2 measurements. This exercise clearly illustrates the influence 
of the occulter size on CME mass measurements and suggests that those measurements, 
where the CME is only partially imaged, in any coronagraph, should be corrected via the 
procedure outlined below. 

Figure [5] shows the 3D CME mass evolution against deprojected height for the 2010 April 3 
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event. The observed mass increase at low heights is due to the gradual appearance of mass 
in the field-of-view (FoV) of the telescope. By comparing observations with theoretical 
considerations of this geometrical effect, we can derive a fit function, which describes our 
observations, and also enables us to estimate the mass hidden behind the occulter. We 
describe the derivation of that fit function, which is based on the following assumptions: 

• The volume of the CME has the shape of a spherical sector. 

• The CME expands self-similary and adiabatically. 

• The CME mass remains constant during the CME propagation from the low corona 
to the first coronagraph FOV. 

We describe the volume of the CME with the formula of a spherical sector: 

rh p9 p2tt 

Vcme = / / r 2 sin9 dr d6 d<p 

Jo Jo Jo (i) 

^h 3 ■ m V ' 

= —(I -cosO), 

where h is the radius of the sphere (height of the CME leading edge measured from the 
center of the solar disk) and 6 is the width of the spherical sector (CME width). Figure [6] 
shows an image of the 2010 April 3 event together with the boundary of a spherical sector 
outlined. Then the occulted CME volume, V occ , is given by Eq. [T] for height h = h occ , the 
effective occultation size: 

V m = ?^(l-cos0). (2) 
Hence the observed CME volume is 

V ohs {h) = Vcme - V occ = y (1 - cos9){h 3 - h 3 occ ). (3) 

Under the assumptions that i) the CME expands self-similarly, i.e., the width 9 remains 
constant during the evolution, and ii) the CME expands adiabatically, the density change 
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can be expressed by 



7 \ 3 

h n 



P(h) = Po ( ^ ) (4) 



m(/i) = p(/i) ^(fc) = m I 1 - ( %^ ) I . (•-») 



with po the density at the beginning of the CME evolution, calculated as y 2 -. The 
parameter m is the initially ejected mass, i.e. 'true' CME mass at the time when the CME 
becomes the first time visible behind the occulting disk. Hence, the initially ejected mass 
mo can be calculated from the visible fraction of mass m(h) which is not occulted. Using 
Eq. [3] and [4] we derive the observed mass at a given height as 

h 

Because the occultation height h occ is a constant, the second term of Eq. [5] becomes small 
with increasing CME height and m(h) converges to m . At heights h = h occ , m(h) = 0, i.e. 
all the CME mass is hidden behind the occulter. At heights below h occ the function m(h) 
of Eq. [5] is not defined. 

A closer look at the later phase of the CME mass evolution (~ 10 to 15 Rq), observable 
only in the COR2 FoV, shows a slight increase in mass, which cannot be explained by a 
pure geometrical effect due to the occulter (Eq. [5]). It could be due to a pile up, i.e. the 
CME sweeps up coronal material as it propagates through the corona. Another possibility 
could be a mass flow from the low corona in addition to mo, when the CME is observable in 
the COR2 FoV. To take into account such effects, we add the term Am(h — h occ ) to Eq. [5] 
to describe a real CME mass increase (in contrast to the geometrical effect). Including this 
term into Eq. [5] we obtain for the CME mass evolution with height 

3 

h 

with Am the real mass increase per height. Whereas the term mo (^f^) 3 in Eq. [^becomes 
smaller with height, the last term increases. At the height h occ , the observed mass is 
equal to zero, since all of the CME mass is occulted. We stress that Eq. [6] describes the 



m(h) = m I 1 - ( ^ ) I + Am(h - h occ ), (6) 
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observed mass as a function of distance but with the knowledge of its fit parameters also 
the 'true' mass evolution can be reconstructed. In Fig. [7] the fit function applied to the 
COR2 observations of our example event is shown as the dotted red line. With this fit 
we find an initial ejected mass mo of 1.68 • 10 15 ± 2.25 • 10 14 g and a real mass increase of 
2.50 • 10 14 ± 2.57 • 10 13 g/R.0. The two fit components (geometrical effect and actual mass 
increase) are plotted separetely. The red solid line shows the 'true' CME mass evolution, 
corrected for the geometrical effect and calculated by 

m(h) = m + Am(h — h occ ). (7) 

For a proper evaluation of the fit function we refer to the appendix of this paper. 



4.1. Mass evolution derived from fit function 

We derive the 'true' mass evolution from a) combined COR1/COR2 measurements 
(applied on 24 out of 25 events - one event has too few data points) and b) solely from 
COR2 measurements (applied on 23 out of 25 events - two events have too few data points). 
Figure [8] shows the results for 10 events of our sample. We combine COR1/COR2 data, 
assuming that for COR2 the occulter effect becomes negligible at larger heights, by using 
all COR1 measurement points and only those COR2 measurements which are at least 10% 
larger than the last measurement in COR1. The combined COR1/COR2 fit curves (Eq. [7]) 
are represented in Fig. [8] as black dotted lines and demonstrate very well that it is possible 
to combine mass measurements derived from different instruments. With this we are able 
to calculate the 'true' mass value from low coronal heights (< 3 R Q ) up to about 20 R Q . 
Comparing the fit parameter for the real mass increase, Am, derived from COR1/COR2 
and COR2 measurements (Eq. [6] shown as red dotted lines in Fig. [8]), we find for about 
50% of the studied events a good match with differences less than 20%. This result is 
expected since we observe in both instruments the same CME. However, for some events 



13 



there are large discrepancies. It turned out that only for well observed events, i.e. if coronal 
disturbances are small compared to the CME observation, the COR1/COR2 fit gives 
reliable results for the real mass increase. Therefore, we determine Am from the fit applied 
to COR2 observations and use the combined COR1/COR2 fit to get an estimate for the 
initial ejected mass mo for which we derive values ranging from 3.6 • 10 14 g to 8.9 • 10 15 
g with a mean value of the logarithmic data of 15.32, corresponding to a mass value of 
2.1 ■ 10 15 g. The distribution of all mo values is displayed in Fig. [9] 

Beside Am, we derived from COR2 the fit parameters m 10 , the mass value at 10 R , and 
m ent i, the mass value at 20 R . For the calculation of m w and m en d we use the 'true' mass 
evolution (Eq. [7]) after eliminating the geometrical effect. The statistical properties of fit 
parameters and CME masses are summarized in Table [T] Since the CME mass values are 



lognormally distributed (see Vourlidas et al. 2010), the statistical quantities are calculated 
in a logarithmic space. 



Figure [10] shows the distribution of the real mass increase rate Am for our sample. Only one 
event out of 23 shows no mass increase in the late phase of the CME evolution. For Am we 
find a range of 2.4 • 10 13 to 2.4 • 10 15 g/R©. The mean value of the logarithmic data is 14.25, 
corresponding to a mass increase of 1.8 • 10 14 g/R©. The relative mass increase Am/m en d 



(see Fig. 11) is on average 0.04, meaning that we have a real CME mass increase of about 
4% per R up to 20 R . Minimum and maximum values are 0.02 and 0.06, respectively. 



4.2. Force and energy distributions 

We use the CME peak accelerations a max , peak velocities v max and the ejected masses 
m to calculate total forces (F = a max m ) and kinetic energies {E^ n = \rriQV^ nax ) in the 
low corona. Since the distributions of the total forces and kinetic energies also follow a 
lognormal distribution, the distributions are plotted and their statistical properties are 
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calculated in a logarithmic scale. The statistical properties are summarized in Table [TJ 
There are different forces (Lorentz force, effective drag force, gravitational force), which 



act on the CME and causing its acceleration or deceleration (Vrsnak 2006). In this work, 
we don't distinguish among the different forces but observe the effect of the net force. 



Figure 12 shows the distribution of the maximal total forces. We note that a max values 
are measured between 1.1 and 3.3 R , whereas mo is assumed as the mass the CME has 
reached in the height range between 1.4 and 2.8 R Q (derived from the fit). But because 



the peak acceleration for many CMEs occurs at low heights (see Bein et aL] 2011 ) that 
they are inaccessible by coronagraphs, mo is the best available estimate for the CME 



masses at that heights. This is confirmed by the Aschwanden et al. (2009) study, who 
did not find large differences in the masses, derived from EUVI and COR1 observations, 



respectively. The distribution in Fig. 12 shows that the derived net force ranges from values 
of 2.2 • 10 14 to 5.1 • 10 16 N. We find for the logarithms of the total forces a mean value of 
15.18, corresponding to 1.5 ■ 10 15 N. 

We also use m to calculate the kinetic energy E kin . v max is measured at heights between 
1.5 and 8.4 R Q , but the velocity does actually not change that much during the evolution 



in the coronagraphic FoVs. In Fig. 13 the distribution of the kinetic energies is plotted, 
which range from 1.0 • 10 23 to 6.6 • 10 24 J. We find a mean value of 23.89 derived from the 
logarithmic data, corresponding to 7.8 ■ 10 23 J. 



5. Summary and Discussion 

In this paper we presented the 3D CME mass evolution against time and deprojected 
heights for a sample of 25 events, which were observed by the COR1 and the COR2 
instrument onboard both STEREO s/c. 

The CME mass evolutions show a strong increase at the beginning of the evolution, which 
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was also found in former studies (Vourlidas et al. 2000 Colaninno & Vourlidas 2009 Carley 



et al. 2012) and is explained by a geometrical effect, i.e. by a gradual appearance of CME 



material above the occulting disk in the course of CME expansion. This occulter effect 
is confirmed by finding a match between temporally overlapping COR1 and COR2 mass 
measurements when applying the same occulter size for both instruments (cf. Fig. [3]). 
However, the observed mass increase in the later phase of the CME evolution (~ 10-15 
R ) is not due to mass hidden behind the occulter but a real increase, which can be 
explained either by pile up at the CME leading edge or continuous mass flow from the lower 
atmosphere. We developed a fit function taking into account both effects of mass increase 
(geometrical and physical) and found a very good match between fitted and observed data 
points. This function returns three fit parameters: (1) the effective occultation size of the 
instrument (depending on the direction and the CME expansion), (2) the initially ejected 
mass of the CME and, (3) the real mass increase per height. In the appendix we show all 
three are physically meaningful parameters. With Eq. [7] we can estimate the 'true' mass 
evolution and thus can derive the 'true' mass value at every height larger than the effective 
occultation size. By this we can calculate m en d, defined as the mass at 20 R and find that 
it is the most appropriate measure for the mass in the late phase of the CME evolution for 
three reasons: i) m end is always related to the same height, ii) m end is derived from the fit, 
which considers the shape of the CME mass-height curve (e.g. a probable mass increase), 
iii) at ~20 R Q the CME reaches the Alfvenic critical point beyond which it can be assumed 



that the evolution more or less ceases (Hundhausen 1972). 

We have shown that mass measurements derived from different instruments can be 
combined even if the coronagraphs have different occulter sizes. We derived from combined 
COR1 and COR2 observations initially ejected masses m in the range of 3.6 — 89 • 10 14 
g. This fit parameter gives us an estimate of CME mass low in the corona (<3 R Q ) which 
may be important for calculating initial forces acting on the CME, since the CME peak 
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accelerations occur also at similar heights low in the corona QBem et al.|2011p . For the total 
forces we found values between 2.2 • 10 14 and 5.1 • 10 16 N. Kinetic energy values were found 
to lie between 1.0 • 10 23 and 6.6 • 10 24 J. We note here that these deprojected quantities are 
in good agreement with the overall CME statistics ( Vourlidas et al. |2010 ) which are based 
on projected quantities (as do most of the relevant work before 2007). This agreement gives 
us confidence that large statistical studies provide robust information for the actual physical 



parameters of CMEs even though they are based on single viewpoint measurements. Emslie 



et al. (2004) showed that CME kinetic energy values can be even higher. They presented 
projected kinetic energy values of 1.8 • 10 25 J and 1.1 • 10 25 J for two strong CME/flare 



events. 



Carley et al. (2012) presented a case study concerning the deprojected energy and 



force evolution of the 2008 December 12 CME event and found values of 6.3 ± 3.7 • 10 22 J 
for the energy and 3.9 ± 5.4 • 10 13 N for the force at about 3 R Q , which is small compared 
to our study. That can be explained by the fact that our sample includes mostly events 
during the ascending phase of the solar cycle and, thus, have on average larger velocities, 
accelerations and masses. 

The derived mass increases per R were found to be between 2% and 6% up to 20 R , 
corresponding to 2.4 • 10 13 g/R© to 2.4 • 10 15 g/R©. It is well known that absolute mass 



measurements are subject to errors of up to a factor of two (Vourlidas et al. 2000 Lugaz 



et al. 2005 Vourlidas et al. 2010). However, here we were mostly interested in the relative 



change of CME masses, not in the absolute calibration, and the relative error is much 
smaller. If the derived mass increases were only due to errors, then we would expect both 



increases and decreases of mass across different events. This is not the case. Lugaz et al. 



(2005) reported from 3D MHD simulations of CMEs a mass increase of the order of 13% 
within 1 to 2.5 hours after the CME initiation and a weaker one up to 1 AU. The authors 
also compared the results from 3D simulations with mass measurements from line-of-sight 
images and found for the latter a weaker mass increase, which they explained by the 
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following: When measuring the mass from a single viewpoint coronagraphic image, it is 
assumed that all mass is placed on the POS, which is not true for the whole volume of the 
CME, leading to an underestimation of the mass. Because of the expansion of the CME, 
this effect becomes larger with height. More electrons deviate from the POS with larger 
distances making it difficult to distinguish a mass increase against the reduced brightness. 
Thus, the mass increases found in our study can be considered as a lower limit. 
These real CME mass increases can be explained by a pile up of mass during the 
propagation/expansion of the CME through the corona, as suggested by simulations from 



Das et al. (2011), who showed that a pile up is possible also at lower heights (2 to 7 R ). 



However, no such effect has been observed so far in the 2.2 - 30 R range (Howard & 



Vourlidas 2005). On the other hand, mass could be supplied to the CME from mass flow 



resulting from the opening of the field and the creation of a temporary coronal hole behind 
the CME in the low solar corona. Coronal dimming regions analysed by former studies 



showed mass outflows lasting for hours (Reinard & Biesecker 2008 Aschwanden et al. 2009 



Bemporad & Mancuso 2010 Miklenic et al. 2011 Tian et al. 2012) and thus it is likely 



that they provide the CME with material also in its later evolution. Vourlidas et al. (2000) 
already pointed out this effect to explain the backward (relative to the CME front) motion 
of the CME center of mass as the event expands in the outer corona. 

The behavior of the center of mass of the CME during its evolution should clarify which of 
the processes described above is responsible for the real CME mass increase observed. In 
case of pile up, more CME material accumulates at the front of the CME and the center 
of mass is expected to move toward the CME leading edge with time. If there is a mass 
flow from behind, the center of mass would drop backward with time. From the COR2 
observations of our events we calculated for each time step the center of mass position 



within the defined CME sector following Vourlidas et al. (2000). The distance of the center 
of mass to the Sun center is then plotted together with the height of the CME leading 
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edge, which we consider as the outer boundary of the modeled CME flux rope (cf. Fig. 14). 
The inner boundary of the flux rope was calculated by the width derived from the forward 



modeling. The distance of both boundaries (black solid lines in Fig. 14 ) increases because 
of the self-similar expansion of the flux rope. Because a significant fraction of the flux rope 
is occulted in the beginning of the CME evolution, the center of mass was not measured 
until the full flux rope could be observed, i.e. the height of its inner boundary is larger than 
the effective occultation size (indicated by the shaded area below the dotted line in Fig. 



14). Fig. 14 indicates with grey lines steps of 10% of the width of the CME body. From this 



we obtain that on average, the center of mass lies on the 25% line and shifts with 5±7% 
towards the inner boundary of the flux rope. For 14 events the center of mass moved with 
more than 5% of the flux rope width toward the inner boundary (mass flow), in 2 cases we 
measured a motion over 5-10% of the CME width toward the outer boundary (pile up). 
For 7 cases the center of mass motion relative to the CME leading edge remains roughly 
constant (i.e. changes are < 5% of the flux rope width) and for 2 events we cannot make a 
clear statement. Therefore, we find that mass flow behind the CME is the main reason for 
the mass increase. 

Prominences, which are often related to CMEs, can also influence the CME mass 
measurements. In some cases their bright Ha radiation can be observed in white light 
coronagraphic observations. When this radiation is wrongly interpreted as Thomson 



scattering, excessive mass values are obtained. Carley et al. (2012) explained the short-time 
increase and decrease of mass in the COR1 measurements of their event under study by 
this phenomenon. Therefore we took special care to exclude events which showed Ha 
radiation in the coronagraphic observations. Nevertheless Ha emission could be one reason 
for the differences in the mass increases derived from the COR2 fit and the combined 
COR1/COR2 fit (see for example 2010 August 27 in Fig. g). Scattering effects of the COR1 
instrument or mass draining down before it reaches the COR2 FoV can also lead to larger 
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mass measurements in the C0R1 FoV. Coronal disturbances, e.g. other CMEs, background 
streamers or particles produced by a solar flare and hitting the detector may influence 
COR1 measurements as well as COR2 measurements. 



6. Conclusions 

The observed CME mass evolution in coronagraphic observations is easily explained by 
the 'true' mass evolution and a geometrical effect. We derived a fit function, which enables 
us to consider both components separately. By this 'true' mass values at every height, 
between the effective occultation size and 20 R , can be estimated. For the majority of 
the events we find a significant CME mass increase, which is in the range of 2 to 6% and 
becomes most important over the distance range 10-20 R . We find that most of the mass 
is located in the rear part of the CME body and, in about half of the events, there is strong 
evidence that the derived mass increase is supplied to the back of the CME by enhanced 
flow from the low corona due to temporary coronal holes created by the CME eruption. 



Pile-up of coronal material ahead the CME is rare in agreement with the findings of Howard 



& Vourlidas (2005) results. 
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mimimum 


maximum 


arithmetic mean ± 
standard deviation 


median ± median 
absolute deviation 


log(m [g]) 


14.56 


15.95 


15.32 ± 0.29 


15.35 ± 0.18 


h OC c [Rq] 


2.69 


5.08 


3.92 ± 0.61 


3.82 ± 0.36 


log(Am[g/R e ]) 


13.38 


15.39 


14.25 ± 0.38 


14.24 ± 0.14 


log(m 10 [g]) 


14.89 


16.45 


15.53 ± 0.34 


15.53 ± 0.23 


log(m end [g]) 


15.11 


16.72 


15.72 ± 0.33 


15.70 ± 0.20 


log(F[N}) 


14.34 


16.71 


15.18 ± 0.56 


15.27 ± 0.41 


log(E kin [J}) 


23.00 


24.82 


23.89 ± 0.48 


24.00 ± 0.42 



Table 1: Statistical quantities derived for different CME mass and fit parameters. m is a 
parameter from the fit applied to combined COR1 and COR2 observations. h occ and Am are 
fit parameters of the fit applied to COR2 observations. m 10 and m en( i are derived from the 
same fit, representing the CME mass at 10 and 20 R , respectively. The force is calculated 

2 

by F = modmax, the kinetic energy is calculated by E kin = — v 2 max . 
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Fig. 1.— STEREO-A (ST-A) and STEREO-B (ST-B) observations of the CME that oc- 
curred on 2010 April 3. The top row shows running difference EUVI images, the second 
and third row COR1 mass images and the last two rows COR2 mass images. The temporal 
evolution of the event can be seen from left to right and top to bottom. The identified CME 
leading edge and POS propagation direction are overplotted in black. The sectors used for 
the CME mass calculations are outlined in white. 
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Fig. 2. — Deprojected kinematics of the 2010 April 3 CME. The top panel shows the height 
of the CME leading edge, measured from the solar surface at the location of the CME source 
region, against time. The plot symbols (crosses for EUVI, squares for COR1 and asterisks 
for COR2) represent the measurement points, the red solid line a spline fit to the data. 
Measurement errors (0.03 R for EUVI, 0.125 R for COR1 and 0.3 R for COR2) are 
plotted but are smaller than the plot symbols in many cases. The first and second time 
derivative of the measurement points and spline fit are shown in the middle and bottom 
panel, representing the velocity and acceleration profile of the CME. 
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Fig. 3. — CME mass evolution with time for the 2010 April 3 event. Blue triangles represent 
COR1 measurements, green triangles mass measurements derived from COR1 observations 
with applying a COR2 occulter. The red asterisks represent COR2 mass observations. The 
top panel shows measurements derived from STEREO-A observations, the middle panel from 
STEREO-B observations and at the bottom panel the 3D CME mass evolution is shown. 



- 25 - 



3-Apr-2010 10:00:00.006 



2000 I 



in 
u 

0) 

o 



>- 



■2000 




-2000 2000 
X (arcsecs) 

3-Apr-2010 09:54:00.00:; 



1 .5x10 



1 .0x10 



5.0x10 



in 

o 
cu 
m 
u 

o 



-5.0x10 



1 .0x10 



1 .5x10 




1.0x10 





X (arcsecs} 



1 .0x10 



Fig. 4. — Mass images from the 2010 April 3 event derived from COR1 (top) and COR2 
(bottom) images recorded close in time. The size of the solar disk is represented by a black 
circle, the size of the COR2 occulter as the green circle. The black square in the COR2 
image shows the FOV of the COR1 image plotted in the upper panel. 
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Fig. 5. — 3D CME mass evolution against deprojected height for the CME observed on 
2010 April 3, derived from COR1 observations (blue triangles) and COR2 observations (red 
asterisks). A fit to the data points together with the fit parameters are plotted in red for the 
COR2 observations and in black for the combination of COR1 and COR2 measurements. 
For the latter fit, COR1 measurements and all COR2 measurements, which are at least 10% 
larger than the last COR1 mass measurement, are used. 
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Fig. 6. — Top: CME image with overlaid boundary of a spherical sector (green). Bottom: 
spherical sector, representing the CME volume (green) with height h and spherical sector 
representing the occulted material (purple) with height h occ . 
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Fig. 7. — 3D CME mass evolution for the 2010 April 3 event. The measurements (asterisks) 
together with the fit (dotted line) and the real mass evolution (solid line) are plotted in red. 
Additionally the two fit components are plotted separetely, the geometrical effect in green, 
the 'real' mass increase as the blue solid line. The black dashed line represents the initial 
ejected mass m . 
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Fig. 8. — 3D CME mass evolution against deprojected heights for 10 events out of our 
sample. Blue triangles represent C0R1 observations, red asterisks C0R2 measurements. 
The fit applied to C0R2 observations is shown as the red dotted line together with its three 
fit parameters. The black dotted line represents the fit applied to combined COR1 and 
COR2 measurements (only COR2 data points were used, which are at least 10% higher than 
the last COR1 data point). From the parameters of this fit (printed in black) the 'true' CME 
mass evolution can be calculated (black straight line). 
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Fig. 9. — Ejected mass m derived from the fit applied to combined COR1 and COR2 
observations. 



- 31 - 




13 



14 15 
log(Am [g/R Sun ]) 



16 



Fig. 10. — Distribution of the real CME mass increase rate Am, derived from the fit applied 
to the COR2 mass measurements. 
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Fig. 11. — Distribution of the relative mass increase per R defined as Am/m end * 100. 
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Fig. 12. — Distribution of total forces calculated by m a max for 24 events. m is derived 
from the fit applied to combined COR1 and COR2 measurements, a max is the peak value of 
the CME acceleration profile. 
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Fig. 13. — Distribution of the CME kinetic energy. 
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Fig. 14. — Center of mass evolution within the CME flux rope for example events. The 
center of mass evolution against the leading edge height is represented by crosses, the outer 
and inner boundary of the CME flux rope are plotted as black solid lines. The geometric 
flux rope center (dashed line) and grey lines in 10% steps of the flux rope width are plotted 
to better track the center of mass motion (crosses) within the flux rope. The dotted line 
reDresents the effective occultation height h„„ derived from the fit. 
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A. Evaluation of the fit 

In the following we evaluate our fit function against the observational parameters. 
The fit parameters are the effective occultation size of the instrument h occ , the initially 
ejected mass of the CME mo and the real mass increase per height Am. As observational 
parameters we use the actual size of the occulter disk and m Q b s . 

Due to the evolving CME mass with height, we define m D ft s as the mean value of the last 
three COR2 measurements, which are measured in the outer range of the COR2 FoV. 



Figure [15] shows the distribution of all derived m t, s values, plotted on a logarithmic scale. 
From the logarithmic data, log(m b s ), we find a mean value of 15.55, which corresponds to 
a mass of 3.6 ■ 10 15 g. Minimum and maximum values of m b s measurements are 1.0 ■ 10 15 g 
and 4.2 • 10 16 g, respectively. 

In the course of the derivation of the fit function we assume that m is the initally ejected 
mass, i.e. a physical meaningful measure. To confirm our assumption we compare m with 



the measured mass values m b s (see Fig. 16). Because m represents the mass at about 



1.5 to 3 R and m b s the mass at about 15 R these two measurements should be highly 
correlated since we observe a linear real mass increase Am. Indeed we find a positive 
correlation with a correlation coefficient of 0.83. 

We also find a positive correlation between m Q b s and the real mass increase Am of c = 0.8£ 



shown in Fig. 17, which means that more massive CMEs also have a larger real mass 



increase. On the other hand, the relative mass increase Am/m en d (see Fig. 11) shows no 
correlation with m b s - 

We derive the effective occultation size h occ from the fit applied to deprojected data, thus 
h occ is larger than the physical radius of the occulter dependent on the deviation of the 
CME propagation direction from the POS. Considering the CME propagation direction, 



h occ gives us meaningful values for the effective size of occultation. Fig. 18 shows the 
distribution of the occultation size projected back on the spacecraft POS. 
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Fig. 15. — Distribution of m j, s , the mean of the last three CME mass measurements derived 
from COR2 observations for a sample of 25 events. 
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Fig. 16. — Ejected CME mass m against m b s , the mean value of the last 3 mass measure- 
ments of C0R2 observations. The regression line is plotted in red. mo is derived from the 
fit applied to combined COR1 and COR2 measurements. 
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Fig. 17. — Real mass increase rate Am (derived from the fit applied to the COR2 mass 
measurements) against m b s , the mean value of the last three C0R2 measurements. The 
regression line (red solid line) and the correlation coefficient c are calculated in double loga- 
rithmic space. 
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Fig. 18. — The effective occultation size h occ is derived from the fit, which is applied to 3D 
COR2 CME masses against deprojected heights. The distribution shows the projection of 
the h occ values on the POS corresponding to the CME propagation direction. STEREO-A 
occultation sizes are plotted in color, STEREO-B occultation sizes in grey. 
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